Isolated Minkowski vacua, and stability analysis 
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We study a recently proposed model, where a codimension one brane is wrapped around the 
axis of symmetry of an internal two dimensional space compactified by a flux. This construction is 
free from the problems which plague delta-like, codimension two branes, where only tension can be 
present. In contrast, arbitrary fields can be localized on this extended brane, and their gravitational 
interaction is standard 4d gravity at large distance. In the first part of this note, we study the de 
Sitter (dS) vacua of the model. The landscape of these vacua is characterized by discrete points 
labeled by two integer numbers, related to the flux responsible for the compactification and to 
the current of a brane field. A Minkowski external space emerges only for a special ratio between 
these two integers, and it is therefore (topologically) isolated from the nearby dS solutions. In the 
second part, we show that the Minkowski vacua are stable under the most generic axially-symmetric 
perturbations (we argue that this is sufficient to ensure the overall stability). 



I. INTRODUCTION 

Models with two extra dimensions have been the ob- 
ject of many studies. Among them, there are the first 
examples of supersymmetric compactifications to four di- 
mensional Minkowski space where a stable spherical 
internal space is achieved through the flux of a gauge field 
and a cosmological constant. More recently, six dimen- 
sional models have been reconsidered in the brane-world 
scenario Q , since two is the minimum number of flat ex- 
tra dimensions to have a fundamental scale of order TeV 
and a compactification of (sub)millimeter range. The 
background solution of [l| can be easily extended to a 6d 
braneworld. If one imposes that the internal and the ex- 
ternal coordinates are factorized, then the presence of a 
brane modifies the internal geometry by creating a deficit 
angle [|[ , in the same way as a point mass modifies a two 
dimensional space In this way, the spherical internal 
space of [l[ is modified to a (so called) rugby ball, with 
conical singularities at the two poles generated by the 
codimension two brane and its Zi image. It is manifest 
from the solution that the only quantity related to the 
brane tension is the deficit angle in the bulk, so that one 
may hope that this construction can provide a self tuning 
mechanism responsible for the vanishing of the 4d cosmo- 
logical constant. Unfortunately, this is not the case for a 
number of reasons [5j]. Nonetheless, the construction of 
[H remains of great interest, since it is an extremely sim- 
ple but complete solution of unwarped extra dimensions, 
which can be employed in a number of applications of 
physical relevance p. 

There is, however, a further problem that these ap- 
plications face, related to the localization of matter and 
gauge fields on the codimension two branes. In field the- 
oretical studies, matter and gauge fields on the brane are 
usually treated as test fields, which do not contribute to 
the background geometry (namely, they are neglected in 
the Einstein equations of the system). However, a more 



complete approach, with gravity also taken into account, 
shows that only tension can be present on the brane, 
while fields with a different equation of state necessary 
lead to worse than conical singularities at the brane loca- 
tion Q. Contrary to what happens for conical singular- 
ities (where the scalar curvature is a delta function sup- 
ported at the tip of the cone) such singularities are not 
integrable, so that the codimension two brane cannot be 
consistently treated even in a distributional sense. This 
appears as a specific example of the problem of defining 
sources of codimension two and higher in general rela- 
tivity H (see also for a recent discussion focused on 
exact codimension 2 solutions). 

A way out of this problem was proposed in 
through the addition of Gauss-Bonnet terms in the bulk, 
and in [TlT |. where the codimension two brane emerges 
at the intersection between two codimension one branes. 
A possibly more conservative approach is to replace the 
delta-like strict codimension two brane with an extended 
codimension one defect [12] (extended defects on a codi- 
mension two bulk were also discussed in [H[). In this 
model, a region close to the singularity is replaced by a 
spherical cap, and a codimension one defect is placed at 
the junction between the two spaces (the same is done 
for both singularities, so to preserve the Zi symmetry of 
the system; moreover, the defect wraps around the main 
axis of symmetry of the system; the overall configuration 
can be seen in fig. 1 of [l2|). Fields with arbitrary equa- 
tion of state can be localized on the defect, and their zero 
modes (that is, the axially-symmetric ones) can be iden- 
tified with the 4d fields of the observable sector. It was 
shown in [l2j that the gravitational interaction between 
brane fields is described by Einstein 4d gravity at large 
distance. Therefore, this construction appears as a phe- 
nomenologically viable and complete theory of two extra 
dimensions, without the need of invoking higher order 
gravity terms, or a more complicated system of branes. 

A natural extension of this construction is the study 
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of more general background solutions. For instance, 
ref. [3] embedded the codimension one brane in a 
Minkowski compactification characterized by a warped 
internal space. Alternatively, one can study cases of cos- 
mological relevance, where the brane is embedded in a 
time dependent background. We perform the first step 
in this direction by studying and classifying the de Sit- 
ter (dS) vacua of this model. More precisely, we look 
for solutions characterized by a dS external geometry, 
and a static internal space. More general solutions can 
presumably be obtained analytically for small deviations 
from the dS or Minkowski ones; in general, the values 
of the bulk cosmological constants and of the brane ten- 
sion determine both the dS expansion rate, and the bulk 
parameters, such as the compactification radius and the 
deficit angle; if more general sources are present, we ex- 
pect a nearly static internal space, and an approximate 
Friedmann- Robertson- Walker cosmology, as long as their 
energy densities are smaller than the brane and bulk cos- 
mological constants (this is typical of extra dimensional 
models with a stabilized internal space). At higher en- 
ergies, we expect a strong evolution of the entire space, 
which can be presumably studied only through numerical 
computations (l5j . 

We find two main results: first, the dS rate must be 
smaller than a given value, related to the inverse of the 
compactification radius. This is in agreement with the 
findings of [l6j |. where it was observed that any com- 
pactification mechanism can stabilize the internal space 
only up to some given expansion rate (by causality rea- 
sons, one cannot expect compactification when the hori- 
zon size becomes parametrically smaller than the size of 
the internal space). A second, more surprising, result, is 
that the Minkowski solutions, obtained for a special re- 
lation between the bulk cosmological constants and the 
brane tension, are separated from the dS ones by topol- 
ogy. Indeed, the possible vacua of the model are labeled 
by two integer numbers. The first integer N is related 
to the winding number of bulk fields charged under the 
gauge symmetry responsible for the compactification of 
the space (if present, these fields impose a quantization 
condition analogous to the one taking place for the Dirac 
monopole). This can also be seen as quantization of the 
flux of the gauge field in the compact space. The sec- 
ond integer n is related to the current of a brane field; 
this current is necessary for matching the discontinuity 
of the gauge field across the two sides of the brane, and 
it controls the position of the brane in the internal space. 
Ref. [13] studied only the Minkowski solutions for the 
system, showing that they occur for N/n = —2. We 
find that a nonvanishing dS rate is instead achieved for 
greater ratios. In the Conclusions, we comment on the 
implications that this can have for the cosmological con- 
stant problem. 

Clearly, for the construction of [l2j to be of any in- 
terest, one should show that it is stable. Only the zero 
modes of the perturbations of this geometry were studied 
in [l2j ■ In the second part of this work, we study the mas- 



sive modes, to ensure that the system has no tachyonic 
instability (we do so only for the Minkowski vacua; from 
what we already mentioned, and from the stability s tudy 
for the spherical compactification without branes (17I ] . 
we expect that the system becomes unstable at high H ) . 
More precisely, we concentrate only on axially symmetric 
perturbations (around the axis of symmetry of the back- 
ground). The reason is mostly technical, since, as we 
will see, already this system of modes is quite involved. 
However, we have a second (although non rigorous) justi- 
fication for this restricted study. The construction is ob- 
tained by cutting in an axially symmetric way the rugby 
ball and the spherical compactifications, and by joining 
them across the brane. Before cutting them, both these 
configurations are stable so we expect that an insta- 
bility - if any - will be related to their interface. For 
instance, the rugby ball may prevail over the spherical 
cap, so that the brane in between them would shrink to- 
wards the pole. Alternatively, the spherical cap may be 
favored, and the string would then extend towards the 
equator. Such instability would show up as a tachyonic 
axially symmetric mode; recently, the stability of 6D chi- 
ral gauged supergravity, including the unwarped "rugby- 
ball" solution, was studied in [la ]. Also that analysis is 
restricted to axially symmetric perturbations; it is ar- 
gued that the study of more general modes is unnec- 
essary, since any angular dependence would contribute 
positively to the corresponding Kaluza-Klein mass. 

We decompose the perturbations into scalar, vector, 
and tensor modes (where the names refer to how they 
transform with respect to transformations along the non- 
compact coordinates) . These three sectors are decoupled 
at the linearized level, and they can be studied sepa- 
rately. The equations for the zero modes studied in [12j 
could be solved analytically. Unfortunately, for massive 
modes, only the tensor mode can be obtained analyti- 
cally, while the bulk equations for the vector and scalar 
modes have to be solved numerically. 1 We do so with 
a shooting method. No tachyonic solution emerged from 
the (rather extensive) numerical computation, so that we 
can conclude 

The paper is organized as follows. In Section 2 we re- 
view the model introduced in [l2j]. Section 3 is devoted 
to the study of the dS backgrounds. In Section 4 we show 
that the system is stable under the the most general set 
of axially symmetric perturbations. The concluding Sec- 
tion 5 contains some remarks on the phenomenological 



1 We are aware of two studies which are close to the present one, 
where the equations in the bulk could be decoupled, and then 
studied analytically. For a spherical bulk, the modes can be 
decomposed on spherical harmonics, and then decoupled due to 
orthogonality in the bulk [l7ll . This cannot be done for the "com- 
posite" bulk solution that we are investigating. Second, in the 
study of |18( the modes could be decoupled by using the equation 
for the dilaton present in the supergravity action. This mode, 
and the corresponding equation, is absent in our case (for other 
studies of linearized gravity in 6d contexts, see Il9l0 . 
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aspects of the model, with a particular focus on the possi- 
ble implications of the separation between the Minkowski 
and dS solutions. 



II. MINKOWSKI COMPACTIFICATION 



We summarize the construction of [121 ], 
the model is 



The action of 
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(1) 
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where Sb is the action of the branc, while S ^ is the 
action in two bulk regions separated by the brane. The 
line element in the two regions is 



ds\ = •q^dx^dx 1 ' 

ds 2 = ^ v dx^dx v + R 2 (3 2 d6 2 + R 2 (3 2 



- R 2 d9 2 + R 2 (3 2 cos 2 6d<t> 2 



"out" 



cos' 



"in" 



(2) 



where here and in the following x denotes the noncom- 
pact directions (notice the choice of the Minkowski met- 
ric; in the next Section we will instead consider a dS 
external space). The brane is at the background posi- 
tion 9 . The region at < 9 < 9, denoted as the "out" 
bulk, is a portion of the so call rugby ball compactifica- 
tion, characterized by the deficit angle 1 — f3 . The region 
9 < 9 < it/2, denoted as "in" bulk, is a spherical cap, 
with the pole at 9 = ir/2. A Z% symmetry extends this 
geometry to the region — tt/2 < 9 < . In addition, the 
background is axially symmetric (around the axis con- 
necting the poles at 9 — ±ir/2). 

The internal compactification is achieved through the 
gauge field configuration 



F 64> = deA^ = M 2 Rf3 cos 9 
and the cosmological constants 



R(3 



R 



(3) 



(4) 



If we suppress the external dimensions, the brane can 
be viewed as a string wrapped around this axis of sym- 
metry. The brane field a acts as a Goldstone boson (v is 
most easily interpreted as the vacuum expectation value 
of a field which breaks the U(l) symmetry on the brane). 
It generates a current which is necessary to provide the 
discontinuity of the magnetic field between the two bulk 
regions. From its own equation of motion, a = ncj), where 
(due to the periodicity of the coordinate) n is an inte- 
ger. The brane position is then found to be (l2| 



arctan 



1-/3 
Rf3q 2 



e v 



(5) 



while the deficit angle in the out bulk is 



l-f3 



2M 4 irsm9 



(6) 



where T is the four dimensional (i.e., after an integration 
along cf>) energy density of the brane. This result, in the 
limit 9 — > 7r/2 (when the brane shrinks to the north pole), 
reproduces the known relation between the tension of a 
codimension two brane, and the deficit angle generated 
by it. 

To conclude the description of the background, we note 
that, if some field, with charge e under the U{1) symme- 
try is present, the deficit angle must satisfy a quantiza- 
tion condition 



(3 



N 



2eM 2 R 



N integer 



(7) 



Such quantization is also known as flux quantization, 
since it can be recast in the form 



$ B = / d9 



2tt 

>F 6c b = —N 
e 



(8) 



where the integral is performed over the entire internal 
space. 

From the Einstein equations, one then finds [l2j that 
this integer must be related to the winding n of the 
brane field a by N/n = — 2 . As we show in the next 
Section, this relation is actually due to the assumption of 
Minkowski noncompact space. Different ratios between 
these two integers result in a dS external geometry. 



III. DE-SITTER COMPACTIFICATION 

We now generalize the Minkowski solution described 
above to the case of a dS noncompact space, character- 
ized by the expansion rate H (namely, we replace r]^ 
by the dS metric in the line elements |2j)). The bulk 
compactification is achieved for 



An 



M 
2R 2 



(l + 9H 2 R 2 ) 



F g< k = M 2 R(3^l - m 2 R 2 cos( 



(9) 



in the out bulk, and 



M A 



2R 2 f3> 



(1 + %(3 2 H 2 R 2 ^ 



F e<j , = M 2 Ri3^\-m 2 R 2 fi 2 cos 9 



(10) 



in the in bulk. We observe that the expansion rate cannot 
exceed the value of 1/ {\fiR\. This is not surprising in 
the light of the findings of [lg|, where it was shown that, 
for any dS compactification, the expansion of the exter- 
nal coordinates has the generic effect of destabilizing the 
internal space. 

We still look for an axially symmetric solutions, so that 
Ag = , while A^ depends only on 9 . Moreover, A^ = 
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at the poles (from regularity), and is continuous 
across the two branes (so that the brane action is well 
defined). This determines in the bulk; the solution 
cannot be provided on a unique chart. In presence of a 
bulk field charged under this U{1) symmetry, a consistent 
solution is possible only when the quantization condition 

2eM 2 R(i\{l -3H 2 R 2 ) 1/2 sin6 

+ (l-3H 2 l3 2 R 2 ) 1/2 (l-smd)}=N (11) 

with N integer, holds (this computation can be per- 
formed exactly as in the Minkowski case; see [IJJ for 
details). Also in this case, the quantization condition 
can be recast in the form ([8|). 

Let us now discuss the brane equations. As for the 
Minkowski case, a — ncj) , where n is an integer. By 
construction, the transverse metric components (g^u and 
g^tf)) are already continuous across the brane. We are 
then left with three nontrivial brane equations (two sec- 
ond Israel conditions, plus Ampere law, relating the dis- 
continuity of Ferf, in the bulk to the current on the brane) . 
With some algebra, they can be recast in the form 



2n 
~N 



(12) 
(13) 



sin 9 ( V 1 -3H 2 R 2 -T) 



yfl - 3H 2 R 2 f3 2 (l - sin 9) + y/l - SH 2 R 2 sin0 

(14) 



1-0 , s 2A * 
tan = — j 

R(3 M 4 



where q = e v , and we have defined, for shortness, 

1-/9 



T = 



(1 - 3H 2 (3 2 R 2 ) 1/2 - (1 - 3H 2 R 2 ) 1/2 (3 



(15) 



To solve the above system of equations, we eliminate 9 
from (|12p and (|14p . and we combine the resulting equa- 
tion with the two bulk expressions for Ao and Aj. We 
obtain 



l-(3 2 



R 2 



l&\q 2 (A, - A ) 



3 (A, - A )' + 2\q 2 (5A s: - 3A ) + 3AV 
IQXq 2 

3 (A, - A ) 2 - 2\q 2 (3A, - 5A ) + 3AV 



H 



- (A, - A ) 2 + 2\q 2 (A, + A ) - X 2 q 4 



where we have rescaled 



2Ao 

M 4 



Ao , 



2Ai 
M 4 



48A<? 2 



A, 



2K 
M 4 



X 



(16) 



(17) 



These solutions are valid only for A < (Aj — A ) jq 2 
(this is because they are obtained by squaring some of 
the above equations). The maximal allowed value leads 



to H = 1/ (V3-R) , which, as we saw from eq. ©, is the 
highest possible value that the system can have for the 
dS expansion rate. 

We see that the Minkowski compactification requires 
the tuning \/A = (\/Aj — VAT) /<Z • For a small deviation 



VA = 



8 



(18) 



eqs. (jl6|) give 

Ao (V^i + \Ao) q 



3/2 



2K 



5 + (5 2 ) 



Ao 2A 3 /2 (VA l ~V^) 



s + o(s 2 ) 



H 2 = 



6 {VM 
+0 (<5 3 ) 



^V7T iq 5 q 2 (Aj + VA~VA^ + Ap) j2 



12 UK, 



(19) 



From either of (fl"2"|) and HU), the brane position then 
satisfies 



tan# 



/A~- \7A~o 5 



+ J-q S + °^ (20) 



The possible values for the above parameters are con- 
straint by the two integer values N and n . Eqs. (jTTJ) 
and ([T3j) give 



N = 2eM7VAi + (i? 



2n 
TV 



3(VA7 



2AoVAl 



-.IP 



/A" 



/AT 



(21) 
0{H 4 ) 



where the expansion in S has been replaced by an expan- 
sion m H 2 through the last of (fT9|) . 

From the first quantization condition, we see that we 
cannot vary A alone without also varying the bulk cos- 
mological constants. However, the Minkowski compacti- 
fication (5 = 0) can be still deformed continuously into a 
dS one. However, once also the second condition is taken 
into account, the Minkowski solution appears to be de- 
tached from the dS ones. Indeed, the choice 2n/N = — 1 
is only compatible with H = , with the only exception 
of the trivial case of Aj = Ao (in which case, 0=1, and 
the brane is actually absent). 2 

We can gain further insight by estimating the param- 
eters entering in eq. (|2"Tj) . Neglecting the subleading H 2 
terms, and for non hierarchical values of the deficit angle 
(that is, and 1 — of order one), both -^AT and \/A~o, as 
well as their difference, are of order 1/R 2 . In addition, 



2 This conclusion actually holds for arbitrary values of H , and not 
just at small <5, as can be seen by studying eq. I)13|l for 2n/N = 
-1. 
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q 2 R does not exceed one (as can be seen in eq. (|20[) - this 
value controls the ratio between the radius of the brane 
and that of the internal space). Therefore, we find 



1 



2n 
TV 



O (R 2 H 2 



(22) 



From the same reasoning, we also see that the expansion 
at small 5 is actually an expansion for R H <C 1 • 

In the concluding Section we comment on the implica- 
tion of these findings for the cosmological constant prob- 
lem. 



includes the Gaussian normal coordinate choice at the 
brane location, which is the most convenient one to in- 
terpret the gravitational effects measured by brane ob- 
servers). These choices do not fix the gauge completely 
(see [12| for details); however, we can still have a general 
(and unambiguous) study if we perform our computation 
in terms of the combination of modes which are invariant 
under the residual gauge freedom. 

The tensor mode h^ u , and the vector modes T^, a M , 
and VJj , are already invariant. For the scalar sector, the 
invariant combinations are instead 



IV. STABILITY OF THE MINKOWSKI 
COMPACTIFICATION 

The goal of this Section is to obtain the massive pertur- 
bations of the system, to verify whether the background 
solution described in Section 2 has tachyonic instabilities. 
For the reasons mentioned in the Introduction, we focus 
on axially symmetric perturbations. 3 The most general 
perturbations of the geometry of this type are 

ds 2 e = (1 + 2$) dl 2 + 2Adldcf> + (1 + 2C) cos 2 Odft + 
+2 (Tp + d^T) dldx^ + 2 (V^ + d^V) d^dx^ + 
{r}^ (1 + 2*) + 2E^ V + E M + V} dx* dx " 

(23) 

where we have defined the "dimensionful angular coordi- 
nates" 



_ / Rf3d6 "in" 
" < RdO "out" ' 



dej} = Rf3dcf> , 



(24) 



and E, 



dvE^ + d^E,,. The vector modes E^, T M , V M 



are transverse, and the tensor mode /i M „ is transverse and 
traceless. The remaining modes are scalar (the denom- 
ination refers to how these modes transform under 4d 
coordinate transformations). Simultaneously, one needs 
to consider the perturbations of the gauge field, 



5A, 



ai, 



8A, 



ai, 



(25) 



where a M is a transverse vector mode. Following the dis- 
cussion of [lH, we fix part of gauge freedom by setting 
Eft = T = V = ai — a' = (here and in the following, 
prime denotes derivative with respect to the rescaled co- 
ordinate I ) . We further impose 



^branc 



$ (0) = o 



(26) 



that is, we require that the brane lies at the unperturbed 
background position, and that gu — 1 there (this choice 



3 Although involved, it is not hard to extend this analysis to gen- 
eral modes. Since the background is axially symmetric, the gen- 
eral dependence of the modes on the angular coordinate can only 
be of the form exp (in</>) , with n integer. 



$ = $ + E" 

C = C -O'tanOE' 



a? + M 2 6' cos BE' 



(27) 



Tensor, vector and scalar modes arc decoupled at the 
linearized level, so we can study the three sectors sep- 
arately. We do so in the next three Subsections. The 
relevant equations were obtained in (T^ |, where the zero 
modes of the system were then studied. The derivation 
of these equations is not repeated here. 



A. Tensor Modes 

The axially symmetric tensor perturbation can be de- 
composed as 



hftu {x, 9) 



h n (0) C^n (x) 



(28) 



where C Kn are 4d Kaluza Klein (KK) tensor modes, 
and h n their wavefunctions in the bulk. Our goal is to 
find the allowed perturbations, and their 4d masses m n . 
The bulk equation 



<9 2 V + h" - e'tanOh' = 



(29) 



(where d 2 denotes the d'Alambertian operator in 4d) can 
be separated in 



d h n dh n j 



(30) 



where we have moved back to the coordinate, and where 
the parameter fi n is different in the two bulk regions: 



m 2 n R 2 /3 2 "in" 



2 d2 J) 



mi R 



out" 



(31) 



From now on, we suppress the index n for brevity, under- 
standing that we are studying one KK mode at a time. 

The bulk equations must be supplemented by a set of 
boundary and parity conditions. First, we require reg- 
ularity at the two poles, imposing that the first deriva- 
tive of h vanishes there (we can impose this condition 
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FIG. 1: i{v(m 2 )) for the specific choice j3 = 0.9 and 9 = 85°. 



FIG. 2: Smallest masses for the tensor modes, for the specific 
choice f3 = 0.9 and for different values of the brane position 



either on the derivative with respect to 9 or I, since the 
two variables are simply related by a constant rescaling) . 
Second, parity considerations impose that the modes are 
even across the equator, h (—9) = h (9) . Finally, we must 
satisfy the junction conditions across the brane, 



hi n (0) = h out {9) 
Kn(0) = KuM or 



, " lm (0) = P^ i (0) (32) 



no 



The bulk equations are solved by the Legendre func- 
tions P u (x) and Q u (x), where x = sin# and v = —1/2 + 
•\/l/4 + n 2 (we denote by v\ and v the values of this pa- 
rameter in the in and out bulk, respectively). The bulk 
solution which is regular at the poles, even across the 
equator, and satisfies the first of ([32]) is 



Pn (M) 



h ut = A 



cos {^y) Pu {x) - ^-sin (jy) ^"oi x ) 



(33) 



where 
A = 



cos 



& tin ft 



Qu (x) 



(34) 

The only undetermined parameter is the mass square 
of the mode, which enters in the two parameters vi^ . 4 It 
can be found by imposing the only remaining condition 
to be satisfied, namely the second of (T3"2"l) . Specifically, 
for any fixed values of [3 and 9, we (numerically) look for 
the roots of 



(35) 



As an example, Fig .fl} shows the behavior of / (m 2 ) 
for the specific choices of (3 = 0.9 and 6 — 85°. As 
can be seen in the figure, there are no tachyonic modes 
in the spectrum (this is the case also for more negative 
values of m 2 than those plotted). We verified that no 
tachyonic modes appear also for several other values of 
the brane position (those reported in figure ©) and the 
deficit angle (We varied (3 = 0.2, 0.9 in steps of 0.1 
). We also notice the presence of the zero mode already 
studied in 

Fig. @ shows instead the mass spectrum for different 
brane positions, and for the specific choice (3 = 0.9. The 
points at 9 = 90° have been obtained for a codimension 
two brane located at the pole. We see that the mass spec- 
trum for the codimension- 1 model converges continuously 
to the one of codimension-2 (we will see that this is not 
the case in the scalar sector). This can be proven analyt- 
ically, from the study of the equations which determine 
the allowed modes. In the limit of 6 — > it/2, the in part 
of the bulk shrinks to zero, so that all the bulk geometry 
is described by the rugby ball, as it is the case for the 
codimension two case; so, the bulk equations converge to 
that of codimension two. Moreover, the boundary/parity 
conditions that we have discussed above reduce to 



dh 



-, =0 



h{-9) = h{9) , as 9 -> tt/2 (36) 



which coincide with those of the codimension two case. 

In the codimension two case, the spectrum does not 
depend on /3 (this is strictly true for axially symmetric 
perturbations); this also emerges from our numerical re- 
sults (not shown here): the dependence of the spectra 
on the deficit angle is very weak for any finite 9 , and it 
disappears as the codimension one brane is shrunk to the 
pole. 



B. Vector Modes 



4 There should also be an overall normalization, which cannot be 
determined by the linearized system; we have fixed it by setting 
h = 1 at the two poles. 



The linearized Einstein and Maxwell equations for the 
vector modes have been derived in [l2j]. The one for T M 
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simply reads 



d Tfj, = 



(37) 



which immediately indicates that this perturbation has 
not massive modes. We decompose the two remaining 
modes as we did for the tensor case, 

a M = ^2 & » > n Qn ^ ' ^ = ^2 w » >™ ( x ) Wn W l M ' 

n n 

Omitting the index n for brevity reasons, the bulk wave- 
functions satisfy the following system of equations 



cfa 
i6 

d 2 W 



da 



1 dW 



= 



,m + f * + = * (39) 



where fx is defined as in eq. (|3ip. 

Due to the parity choice of the background, the mode 
a must be odd across the equator, while the mode W 
must be even, 

a(-9) = -a(6) , W (—9) = W (9) (40) 

In addition, there are regularity conditions at the poles, 

da _ dW 
d9 ±tt/2 dO |±tt/2 

and junction conditions across the brane. 







(41) 



(0) (9) , W in {6) = W in {0) 

da a \ „ / da a 



d9 tan 9 

[dW „ 

— — + 2a cos ( 
\ a# 







|6 ,in 



Vd6> ta,n 9 



\e ,out 



= /3 



+ 2a cos 



(42) 



\e ,out 



For massless modes, the bulk equations (|39|) form a 
system of two coupled first order differential equations 
in terms of da/d9 and dW/d6. This system can be 
solved analytically. However, for nonvanishing mass, 
these equations must be integrated numerically. 

Therefore, to find the spectrum of vector modes, we 
resort to a shooting method, which is appropriate for 
boundary value problems. Each mode is specified by its 
mass, and by a series of parameters which determine the 
initial conditions at one of the poles. For definitcncss, 
we start from the south pole. As we discuss below, we 
actually need only one such parameter, which we denote 
by C. We start from some guessed values for to 2 and 
C, and we then solve the bulk equations (|39|) (when we 
cross a brane, we impose the conditions ([4"2"| ). If the 
resulting solution turns out to be regular, and to have 
the correct parity assignment across the equator, then 
we have managed to identify one physical mode of the 
system. 

In practice, the bulk solutions that we obtain numer- 
ically never satisfy these properties, signaling that the 



initial guess for the parameters m 2 and C was wrong. 
We can define the two "distances" , 



di 



a(-e)+a(6) \ = / W(-9)-W(9) \ 

a{~9)-a{9)) 2 \W(-9) + W(9)J 

(43) 

which indicate how far the solution is from being Z2 sym- 
metric. We then proceed in two steps: (i) we densely 
scan the parameter space {to 2 , C} within some given 
range; the values leading to the smallest distances are re- 
garded as our best guesses; (ii) we use a Newton's method 
to find the zeros of these distances, starting from the 
best guesses. Provided the initial conditions are dense 
enough, Newton's method converges to all the physical 
solutions of the system, having values of {to 2 , C} not too 
far from the probed range of values. Indeed, the two- 
dimensional nature of the initial parameter space, and 
the fact that the bulk geometry is regular everywhere, 
make the numerical problem a relatively simple one. It 
is easy to verify (for instance, by increasing the density 
of the initial scan) that all the solutions are reached with 
this method. 

The main numerical difficulty occurs at the south pole, 
where the coordinate system used is singular. To over- 
come this, we actually solve (by Taylor expansion) the 
bulk equations analytically in a neighborhood of the 
south pole. As for the tensor sector, there is one overall 
normalization which cannot be determined by these lin- 
earized equations. We fix this by imposing a (— 7r/2) = 1 
. 6 We then find 



a(-7r/2 + e) = 1 + C e 2 + O (e 4 ) 



W(-TT/2 + e) = ( 2C+ r — l < 



0( e 4 ) 



(44) 



The solutions of a system of second order equations, are 
usually specified by the values of the functions and their 
first derivatives at a given point. In the present case, due 
to the coordinate singularity at the poles, we also need to 
specify one of the second derivatives (we also note that 
the linear terms in the expansions vanish, due to the 
regularity conditions (I41|) ). We started our numerical 



More accurately, the vanishing of d\ and d.2 is a necessary, but 
not sufficient, condition for a mode to be a physical perturba- 
tions of the system; indeed, in some cases those conditions were 
(accidentally) satisfied, although the modes did not have the cor- 
rect parity all throughout the bulk. The easiest way to solve this 
problem (in an automated way) is to check the parity condition 
also at other bulk positions. Specifically, we discharged all those 
solutions which satisfied d\ = = 0, but which had the wrong 
parity at the equator (physical solutions satisfy W = da/dO = 
at 8 = 0). In all the cases we attempted, this was enough to 
eliminate all the spurious solutions (we always verified, by direct 
inspection, that all the modes which passed this second check 
had the correct parity all throughout the bulk). 

6 This choice does not include the possibility of a = at the pole. 
In this case, the numerical computation only leads to the zero 
mode characterized by constant a and W in the entire bulk. 
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FIG. 3: Smallest masses for the vector modes, for various 
brane positions and for the specific value (3 = 0.9 . 



C. Scalar Modes 

There are four scalar gauge invariant combinations, 
satisfying the set of linearized equations derived in • 
Two bulk constraint equations (containing at most first 
order derivatives) can be used to express C and in 
terms of the other two modes $ and . This leaves us 
with the two bulk equations 



d 2 $ 


- tan 6 


4™- 




~m ' 


^ de 


3 ^J + 












- tan 6 




d6 2 





2 



2 $ 



2 



f = 



(45) 



evolutions with e = 1 6 , with the values of the wave 
functions and their derivatives obtained from (f44|) . 

We performed the analysis for several values 
of P (namely (3 = 0.3,0.6,0.9) and 6 (8 = 
60,65,70,75,80,85). The initial region scanned was at 
-100 < m 2 < 100 and -20 < C < 20 (the density 
was progressively increased until the final roots did not 
change. The final density we used was such that the /x 2 
values are varied in steps of 0.2 and C values are var- 
ied in steps of 0.02 in every iteration). In our opinion, 
this range is sufficiently large to probe the stability of 
the model against tachyonic modes. Indeed, the first KK 
modes are expected to have a mass of the order of the 
inverse compactification radius, corresponding to \/i 2 \ of 
order one. Also the parameter C is naturally expected 
to be of order one (since it comes from a Taylor expan- 
sion). Moreover, Newton's method can converge to so- 
lutions outside this range of starting values for m 2 and 
C (this is indeed what happened in some cases). In all 
the cases studied, the two distances d\ and d2 strongly in- 
creased at negative m 2 (they were typically several orders 
of magnitudes greater than for positive m ), indicating 
that no tachyons are present in the model; to have a fur- 
ther check, we started Newton's method also from some 
of the guessed values with negative m 2 (despite the cor- 
responding distances d\ and c?2 were always very high); 
the method never converged to any tachyonic mode. 

In figure [3] we show the obtained spectrum for several 
brane positions and for the specific choice (3 = 0.9. We 
note the presence of one massless mode. As for the tensor 
case, we also show the results for the strict codimension 
two case (9 = 90°). We observe that, also for the vec- 
tor sector, the limit of shrinking the extended brane to 
a codimension two defect is continuous (this can also be 
proven analytically from the study of the linearized equa- 
tions, in the same way as we did for the tensor modes). 
Moreover, also for the vector sector, we observed a weak 
dependence (not shown here) of the spectrum on the 
value of the deficit angle (3 (the reason is the same as 
for the tensor sector). 



where fi is related to the physical mass as in ([3"Tj) . 

The parity assignment of the background imposes that 
both modes are even. Moreover, regularity at the poles 
requires 



|±tt/2 



d % _dV _ d$ 



d0 



d6 



Once we insert these conditions in one of the constraint 
bulk equations (which is legitimate, since the involved 
quantities are continuous as we approach the poles), we 
find 



(47) 



$ + * ) \±n/2 = 



The junction conditions at the brane location were also 
expressed in [12| in terms of all 4 gauge invariant scalar 
combinations, plus the the quantity E' . This quantity 
can be identified with the (scalar) perturbation of the 
brane position. In an arbitrary gauge, the brane is at the 
perturbed position — Q + Q (x^) . This quantity changes 
when we performed a change of coordinates involving 
the bulk coordinate I . The combination which does not 
change under such change of coordinate is ( = ( — E' , 
which can be then interpreted as the gauge invariant per- 
turbation of the brane position. Not surprisingly, this 
is the quantity which enters in the junction conditions, 
when they are written in terms of the gauge invariant 
perturbations (p2"T)). We further restricted the gauge free- 
dom by choosing a system of coordinates where the brane 
remains at the background position, that is C, — . In 
this case, C, = —E' , which is the quantity entering in the 
junction conditions. 

We can combine the junctions conditions given in [l2| 
to eliminate E' . This requires the use of the bulk equa- 
tions (which is however legitimate, since the junction 
conditions relate bulk quantities at the two sides of the 
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FIG. 4: Smallest masses for the scalar modes, for various 
brane positions and for the specific choice /3 = 0.9 . 



brane). After some algebra, we find 
[*]./ = 

[(1 + sin 2 0)$ + sin^ cos - $')],/ = 



0' (4*'-/i 2 cos6»sin6l(*-$; 



-5tf'-tan0$ + *') 9'- 

cos sin 9 



(48) 



where [f]j = / ou t — /in denotes the difference of the quan- 
tity / between the two sides of the brane. 

Also for the scalar sector, the bulk equations ([45]) must 
be solved numerically. We therefore perform a numerical 
analysis analogous to the one done for the vector modes. 
We first solve the bulk equations analytically in a neigh- 
borhood of the south pole. Fixing the overall normaliza- 
tion by setting "F = 1 at the south pole 7 , and taking 
into account the regularity conditions mentioned above, 
we find 



*(-7r/2 + e) = l + Ce 2 + 0(e 4 ) 



$(-7r/2 + e) 



AC- 1 



O 



(49) 



Also in this case, the mode is uniquely determined by 
the two parameters m 2 and C. The numerical investiga- 
tion then proceeds in the same way as for vectors. Fig. 0] 
shows the lightest masses in the spectrum, for the specific 
choice of (3 = 0.9 and for different brane positions (the 
small "oscillatory" behavior of the eigenmasses visible in 
the figure is probably due to numerical errors, and it 
gives a measure of the precision of the computation) . As 
for the the other two sectors, the computation does not 



7 This choice does not include the possibility of \E> = at the 
pole (such modes could in principle exist, since they could have 
a nonvanishing second derivative at the pole). We performed a 
separate numerical investigation for this case, which however did 
not show the existence of any such mode. 



show any tachyonic modes (the modes exhibit a very bad 
parity for all negative values of to 2 we have attempted). 
However, there are two interesting differences between 
the scalar and the other two sectors. 

The first difference is the absence of scalar zero modes 
(as can be also verified by solving the equations analyt- 
ically, which is possible for vanishing mass). This indi- 
cates that all the moduli of the model have been lifted 
(by the fluxes and tensions in the system), and that the 
compactification is stable (this is the case also for a codi- 
mension 2 brane at the pole 20]). The present stability 
analysis is done in absence of ( matter or gauge) fields lo- 
calized on the brane. Ref. 12[ studied the gravitational 
interaction between brane sources; it was found that two 
zero modes are then excited, and contribute to reproduce 
Einstein 4d gravity at large distances. A similar situa- 
tion is also encountered in 5d models, for instance the 
Randall-Sundrum model with a single brane [211 ]. The 
background solution of [2l[ has no scalar perturbations; 
however, when (matter) fields are localized on the brane, 
a scalar zero mode - often denoted as brane bending [22j 
- is excited, and gives a relevant contribution to the grav- 
itational interactions between the brane fields. 

The second peculiarity of the scalar sector is that the 
limit 9 — ► 7r/2 is discontinuous. To see this, we solve the 
linearized equations when the brane is close to the pole, 
at the position 9 = —ir/2 + e. The analytical solution 
([49)) then accurately describes the modes in the in bulk 
immediately before the brane. We can then expand the 
junction conditions (148|) for small e, and obtain the values 
of the modes in the out part of the bulk immediately after 
the brane. They are 



d6 

d$ 
d9 



1 + O (e 2 ) 

8C + 3Q3- l)/3m 2 
4^ 

-2 + /3 + 0(e 2 ) 

2 J^l + o { e) 



O (e 3 ) 



(50) 



The limit e — > (i.e. 9 — > 7r/2) would be continuous if, 
these values converged to those which must be imposed 
for a codimension two brane at the pole. The latter values 
are * + l> = d^/d9 = d<P/d9 = , which clearly shows 
that the limit is not continuous (the only exception is 
the trivial case of a vanishing deficit angle, (3 — 1, when 
both cases collapse to a spherical compactification with 
an empty brane at the pole). 

This discontinuity, however, does not lead to any ap- 
preciable discontinuity on the lowest eigenmasses, as can 
be observed from fig. 2) (the values for 9 — ir/2 refers to 
the codimension two brane). For small e , the eigenfunc- 
tions, although starting from a different value on the out- 
side bulk, quickly approach the ones of the codimension 
two case, leading to nearly identical eigenmasses (within 
the accuracy of the numerical computation) . 

We actually observe from the last of ([50)) that the limit 
9 — > tt/2 actually leads to a divergent derivative of 3> on 



10 



the outside bulk. This results in divergent terms in the 
linearized Einstein tensor. An analogous result is also 
encountered for the scalar modes excited by brane fields. 
It was found in [l2| that one scalar mode diverges when 
the codimension one brane is shrunk to the pole. This 
singular limit in the scalar sector is what precludes the 
localization of matter fields on a strict codimension two 
brane. 



V. CONCLUSIONS 

We studied a brane-world model in a six dimensional 
space-time, in which two of the dimensions are compact- 
ified by a flux. The model is characterized by a codi- 
mension one brane, with one dimension extending inside 
the compact space. This construction avoids the singu- 
larities which plague codimension two and higher defects, 
where only tension can be localized. Indeed, it was shown 
in [l2| that fields with arbitrary tension can be localized 
on this defect, and that their gravitational interaction is 
described by Einstein 4d gravity at large distance. 

We do not regard this construction as a regularized 
codimension two brane, in the sense that one does not 
recover a well-behaved solution as the extended defect is 
shrunk to a zero size in the bulk. More accurately, the lin- 
earized computation of gravity in this model breaks down 
(in the scalar sector) as the size of the brane decreases. 
While the only safe claim that one can make is that grav- 
ity becomes strong in this limit, it is hard to expect a reg- 
ular nonperturbative limit, in the light of the fact that 
the strict codimension two case is itself badly singular. 
The positive aspect of this statement is that this con- 
struction leads to distinct and potentially testable pre- 
dictions for short-range gravitational and electroweak in- 
teractions (corrections to standard gravity would show up 
at distances comparable to the compactification radius, 
while electroweak effects would appear at energies close 
to the inverse size of the brane in the internal space). 
We expect that other inequivalent regular constructions 
are possible, but that they would lead to different short 
range observables. 

To make clear predictions, one needs to obtain the mas- 
sive spectrum of perturbations of the model. While only 
the massless modes were studied in [l2j], we performed 
this computation in the second part of this note. Even 
more importantly, the computation is mandatory to ver- 
ify that no tachyonic mode is present, so that the con- 
struction is stable. As also argued in [la ], the study of ax- 
ially symmetric perturbations should be enough for this 
check (since the background is itself axially symmetric, 
modes with a nontrivial angular dependence are expected 
to have a higher mass). For this reason, we focused our 
investigation on these modes. 

Due to technical difficulties that we have discussed in 
the paper, we were not able to decouple the system of 
bulk equations for the vector and scalar modes, so that 
we had to resort to a numerical investigation. We did so 



with a shooting method (slightly modified, to cope with 
the coordinate singularities at the poles). Clearly, nu- 
merical methods can only guarantee the stability within 
the range probed. However, we conducted a rather ex- 
tensive search. While Kaluza-Klein masses are naturally 
expected to be of the order of the inverse compactifica- 
tion radius, we densely investigated a parameter space in 
the interval 8 -100/.R 2 < m 2 < 100/i? 2 , and for several 
bulk parameters (brane position and deficit angle). In 
no case we found evidence for tachyonic solutions. Since 
this is a relatively easy numerical problem (the bulk is 
regular, and there are no strong hierarchies present), we 
believe that the present analysis ensures the stability of 
the construction. 

While the above considerations are valid for a 
Minkowski external space, in the first part of this analy- 
sis we studied the dS solutions of this model. We found 
that the space of vacua is characterized by discrete points 
labeled by two integer numbers, related to the quantized 
values of the flux in the bulk (N) and of the current of 
a brane field (n), which, in turns, controls the position 
of the brane in the internal space. The Minkowksi com- 
pactification requires N = —2n. If this is the case, we 
can actually have different Minkowski compactifications, 
provided the cosmological constant on the brane (A) and 
on the two sides of the bulk (A^o) satisfy 



A, - VAo 



A, 



2e 2 M 8 
N 2 



(51) 



As long as these relations are satisfied, a change in the 
brane and bulk tensions leads to a different internal 
space, but only to Minkowski external geometry. Since 
the ratio of the two integers N and n cannot be varied 
continuously, it is tempting to ask whether this can be of 
some relevance for the cosmological constant problem. 

Discrete vacua typically arise in presence of fluxes, and 
several studies have already attempted to use this as a 
solution of the cosmological constant problem. The orig- 
inal mechanism of [23| is realized with a 4— form in four 
dimensions. This form can acquire only quantized values 
in units of a charge q, and its energy density behaves as a 
cosmological constant. The value of the form can change 
through membrane nucleation; this is however a very slow 
process, and it could be possible that the present universe 
is trapped in a metastable state, where the vacuum en- 
ergy and of the 4— form add up to the observed value of 
the cosmological constant A to t- This mechanism allows 
for several possible values of A to t, and, provided these 
values are densely packed together, one may hope to re- 
produce the observed expansion rate for some value of 
the quantized flux (even if this is not the case initially, 
one should simply wait until the flux tunnels to a value 
compatible with observations). This, however, requires 



8 We actually employed a Newton method which could converge 
to solutions also outside this interval, if present. 
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very small values of q, which - besides being unnatural 
- do not lead to the known cosmology (they would lead 
to a reheating temperature much smaller than the one 
needed for primordial Nucleosynthesis, see [24| for a dis- 
cussion). String theory can offer a big improvement in 
this respect. There, one requires an internal space which 
can be stabilized by fluxes 25], and which can have a 
complex topology, with several non-contractible cycles. 
There are various ways of wrapping fluxes around such 
internal space, leading to several quantized fluxes, and 
to a multi-dimensional set of discrete vacua. Due to the 
high dimensionality, it is natural to expect vacua with a 
value of Atot compatible with the observed one, even if 
the distance between different vacua (namely, the value 
of the charge q) is large [23 |. This is one of t he p ossible 
realizations of the landscape of string theory [26J . 

In the present context, inserting the measured value of 
H in eq. (|22|) . we find 

i+i-w^v (52) 

TV \0.1 mm / 

In the mechanism of [23| . the present value of H is 
achieved at the price of an unnaturally small charge, and 
very large flux. Our realization does not put a significant 
constraint of the charge. However, the value (f52|) requires 
a very tiny "mismatch" of the relation 2n/N = — 1, which 
can be achieved only when the two winding numbers 
are themselves O (lO 60 ). Although the corresponding re- 
quest (very large flux) is usually not listed as a drawback 
of [23] , it is hard to regard such high windings as natural. 

We can think of some possible improvement. The ne- 
cessity of large windings is probably due to the extreme 
simplicity of the model. As we mentioned, the prob- 
lem of localizing fields in general relativity is present for 
any defect of codimension higher than one, and not just 
for codimension two. Assuming that such a construction 
must be done also for more realistic (and richer) mod- 
els, we can expect that the presence of more fluxes can 
allow for a solution without too large winding, in a sim- 



ilar way as the string realization improves over the one 
of 24] . Alternatively, we may be satisfied in providing at 
least a partial solution to the cosmological constant prob- 
lem. Rather than requiring large windings, it is proba- 
bly more natural to assume that we are locked in one 
of the Minkowski vacua, with TV = — 2n. This relation 
(which by itself does not appear unnatural) could pos- 
sibly explain why the "big" cosmological constant van- 
ishes. The coincidence problem could instead be solved 
by some additional field, which would then play the role 
of quintessence. This would still be an improvement with 
respect to the usual models of quintessence, where the 
absence of a "big" cosmological constant is typically left 
unexplained. 

To improve over these considerations requires a better 
understanding of the background solutions of the model. 
For instance, it may be possible to have TV = — 2n, 
even if the cosmological constants do not satisfy the re- 
lations (|5ip . This may be compatible with a more gen- 
eral solution than the restricted ansatz (Minkowski or 
dS external geometry, times a static internal space) as- 
sumed here. Moreover, in order to study quintessence in 
this context, one needs to include sources with a differ- 
ent equation of state than vacuum. It is usually hard, 
if not impossible, to obtain analytical solutions with a 
time evolving internal space. However, such questions 
can be possibly addressed analytically at low energies, or 
numerically along the lines of [151 ]. 
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